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ABSTRACT 

Context. The radius of an exoplanet is one of its most important parameters. Studies of planetary interiors and their evolution require 

1% precision in the radius determination. Transiting exoplanets offer a unique oppurtunity to measure the radius of exoplanets in 

stellar units. These radius measurements and their precision are strongly affected by our knowledge of limb darkening. 

Aims. We study how the precision of the exoplanet radius determination is affected by our present knowledge of limb darkening in 

two cases; when we fix the limb darkening coefficients and when we adjust them. We also investigate the effects of spots in one-colour 

photometry. 

Methods. We study the effect of limb darkening on the planetary radius determination both via analytical expressions and by numerical 
experiments. We also compare some of the existing limb darkening tables. When stellar spots affect the fit, we replace the limb 
darkening coefficients, calculated for the unspotted cases, with effective limb darkening coefficients to describe the effect of the spots. 
Results. There are two important cases. (1) When one fixes the limb darkening values according to some theoretical predictions, 
the inconsistencies of the tables do not allow us to reach accuracy in the planetary radius of better than 1-10% (depending on the 
impact parameter) if the host star's surface effective temperature is higher than 5000K. Below 5000K the radius ratio determination 
may contain even 20% error. (2) When one allows adjustment of the limb darkening coefficients, the a/R, ratio, the planet-to-stellar 
radius ratio, and the impact parameter can be determined with sufficient accuracy (< 1%), if the signal-to-noise ratio is high enough. 
However, the presence of stellar spots and faculae can destroy the agreement between the limb darkening tables and the fitted limb 
darkening coefficients, but this does not affect the precision of the planet radius determination. We also find that it is necessary to fit 
the contamination factor, too. 

Conclusions. We conclude that the present inconsistencies of theoretical stellar limb darkening tables suggests one should not fix 
the limb darkening coefficients. When one allows them to be adjusted, then the planet radius, impact parameter, and the a/Rs can be 
obtained with the required precision. 
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1. Introduction 

The depth of an exoplanetary transit is primarily determined by 
the ratio of the radius of the transiting exoplanet to the host star. 
' Limb darkening is a second-order, but key effect for determin- 
, ing the exact radius of the planet, because limb darkening not 

■ only modifies the shape of the transit light curve, but it also 
significantly affects the true transit depth. There are other third- 
order effects (e.g. night-side radiation, exorings, gravity darken- 

■ ing, etc.), but they are not studied here. 

Transit parameters are determined by fitting a model to the 
observed light curve data points. Typically the number of free 
parameters is around seven (e.g. epoch, period, duration of the 
transit, the impact parameter, radius ratio of the planet and the 
host star, and two limb darkening coefficients). If we knew a 
priori the exact values of the limb darkening coefficients, then 
the number of free parameters would be significantly reduced by 
fixing the limb darkening coefficients. This would mean not only 
that the dimensions of this optimization problem are reduced, 
but also that the occasional degeneracy of the fitting procedure 
reported by several authors would disappear (e.g. Brown et al. 
2001; Deleuil et al. 2008; Csizmadia et al. 201 1). 

Values of the limb darkening coefficients have been pub- 
lished by several authors for many photometric passbands as a 



function of the stellar effective temperature, metallicity, \ogg, 
and turbulent velocity in the stellar atmosphere. Such tables can 
be found in van Hamme (1993), Diaz-Cordoves et al. (1995), 
Claret et al. (1995), Claret (2000, 2004), Claret & Hauschildt 
(2003), Barban (2003), Sing (2010), Howaith (201 1), and Claret 
& Bloemen (2011). For earlier tables see the summary of van 
Hamme (1993). 

Some investigators of transit light curves have found good 
agreement between theoretical predictions and fitted limb dark- 
ening coefficients (e.g. in CoRoT-8, Borde et al 2010, or in 
CoRoT-1 1 Gandolfi et al. 2010); however, sometimes larger dif- 
ferences were found (e.g. at 2cr level in CoRoT-13, Cabrera et al. 
2010) or the fit did not constrain the values of the limb darken- 
ing coefficients at all (e.g. in CoRoT-12 Gillon 2010). Complete 
disagreement can often be found between the theoretically pre- 
dicted and the observed limb darkening coefficients, even where 
the quality of the light curve is extremely good (e.g. the analy- 
sis of the HST Ught curve of HD 209 458 by Claret (2009), re- 
analysis of Kepler-5b by Kipping & Bakos (201 1)); or else, the 
ground-based photometry of WASP-13 by (Barros et al. 2012), 
who also point out that these discrepancies are not caused by the 
photometric signal-to-noise ratio (S/N), and the quality of the 
light curves are sufficient to say that aforementioned discrepan- 
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cies are real. Claret (2009) concludes that uncertainties in the 
stellar parameters (effective temperature, surface gravity (log g), 
metallicity) are also not responsible for the discrepancies. 

Neilson & Lester (201 1) point out that the used plan-parallel 
theory might be replaced by spherical symmetry for calculating 
limb darkening coefficients, which may provide better predic- 
tions for limb darkening. The 3D stellar atmosphere models also 
offer an oppurtunity to improve the predictions. However, these 
very new models have not been tested against full and multi- 
colour light curves yet (Hayek et al. 2012). 

Since limb darkening changes true transit depths (see 
Appendix), the question arises: do these insonsistencies affect 
the accuracy of the transit parameter determination? As we show, 
the answer is yes. In addition, to distinguish between different 
planetary models and to study them in detail we require ±1% 
precision in planetary radii below five Earth masses (Wagner et 
al 2011; F. Sohl, priv. comm.), a level that cannot be reached in 
all temperature regions when fixing the limb darkening coeffi- 
cients for the transit fit. 

The aim of this paper is not only to take a step toward better 
understanding what can cause inconsistent observed limb dark- 
ening coefficients. We also want to make certain that such in- 
consistencies between observed and predicted Umb darkening 
coefficients do not affect the determination of transit parameters. 

In Section 2 we study the differences between recent limb 
darkening tables. In Section 3 we investigate the effect of the 
poor knowledge of the limb darkening on the planet-to-stellar 
radius ratio. In Section 4 we check the accuracy of the adjusted 
limb darkening coefficients and study this on synthetic data. We 
call attention to the impact of stellar spots on the transit fits and 
their importance in the observable limb darkening coefficients 
in Section 5. We check the determinability of the Umb darken- 
ing coefficients in Section 6. Our conclusions can be found in 
Section 7. 



2. Comparison of theoretical predictions for limb 
darkening coefficients 

In Figure 1 the limb darkening coefficients predictions of Sing 
(2010) and Claret & Bloeman (2011) are plotted for three dif- 
ferent stellar metallicities: [M/H] - -0.5,0.0, 1.0. It is obvious 
that the values given by these different authors deviate from each 
other by a value of 0.01 - 0.3, depending on which temperature 
region we consider. The usual range of the limb darkening coef- 
ficients is between approximately -0.1 and 0.8, so these differ- 
ences are considered as large. 

We concentrate on the 3500-8000 K temperature range, 
because most of the transiting exoplanets have been de- 
tected around main sequence stars with such temperatures. 
Interestingly, we have better agreement between models in the 
high-metalUcity range {[M/H] = 1.0). Here the coefficients 
agree with each other, within 0.02 at the low-temperature end, 
but over 7000 K, they differ from each other by 0.05. Fortunately, 
this temperature region is less interesting for transit hunters. At 
intermediate temperatures the predictions for limb darkening co- 
efficients do not differ by more than 0.01 to 0.02. 

At lower metallicities the situation is completely different. 
Below ~ 4500K, the two lower metallicities are in complete dis- 
agreement, and the relative differences between the predictions 
can be as high as 0.3. Between ~ 5000 K and ~ 8000 K the sit- 
uation is still problematic, and in the case of [M/H] = 0.0 the 
two tables generally differ by 0.2. Similar conclusions have also 
been reached by Claret & Bloemen (201 1). 



Some modellers fix the limb darkening coefficients accord- 
ing to these or other tables when they fit a transit curve. As is 
clear from the aforementioned discrepancies, in some tempera- 
ture regions this is not a good strategy because the values of the 
limb darkening coefficients are not consistent and it is not clear 
which are the correct ones. 

When investigators adjust the Umb darkening coefficients, 
then it is customary to compare the results to some limb darken- 
ing coefficient table. Unfortunately, most authors compare their 
results to only one table, and this may not be appropriate. It is 
necessary to compare the results to all available tables because 
of the disagreements between tables. 

Such a comparison was done by Southworth et al. (2007), 
among others. After modelling the light curve of an eclipsing 
binary star, they compared the resulting Umb darkening coeffi- 
cients to the tables of Diaz-Cordoves & Claret (1995), Claret 
(2000), and Claret & Hauschildt (2003). The predicted values of 
these three tables typically differ from each other by ~ 40 - 60% 
(see Table 5 of Southworth 2007). In their paper the relative error 
of the resulting limb darkening coefficients was ~ 35%, and the 
limb darkening values observed were consistent with more than 
one theory. Consequently, they could not decide which theory 
describes the reality, so the uncertainty limits are inconclusive 
for distinguishing the theories. 

Claret (2009) concluds that no presently available Umb dark- 
ening calculation is able to reproduce the surface brightness dis- 
tribution of HD 209 458. Claret (2008) also points out that the 
observations do not support any of the Umb darkening tables. 

The differences between the different Umb darkening tables 
are mainly due to the adopted numerical methods and not to the 
atmosphere models themselves. In addition, to have a clue to 
stellar atmosphere studies, it will be necessary to compare the 
differences between the observations and theory. This gives an- 
other argument for why it is worthy and necessary to fit the Umb 
darkening coefficients. 

3. The effect of poor knowledge of fixed limb 
darkening coefficients on the radius 
determination 

Here we investigate the hidden inaccuracy caused by fixing the 
limb darkening coefficients during the transit-fitting procedure. 
Our goal is to calculate the error in k caused by the theoretical 
uncertainties of the limb darkening coefficients. We assume the 
time being that the normalized transit depth AF/F is perfectly 
known since we are only interested in the error caused by the 
inconsistent limb darkening tables. 

Let Ami and Au2 be the uncertainties of the theoretically cal- 
culated linear and the quadratic limb darkening coefficients, re- 
spectively. These uncertainties have four sources: (i) the theo- 
ries have some uncertainties. These kinds of uncertainties are 
not well known, but as a first approximation, we can assume that 
their order of magnitude is equal to the differences between dif- 
ferent theories as mentioned in Section 2. (ii) The input stellar 
parameters (Tef/, logg, metallicity) have uncertainties that wiU 
be reflected in the prediction of Umb darkening coefficients. The 
uncertainties of stellar parameters usually cause small uncertain- 
ties in the prediction of the limb darkening coefficients, and for 
instance, Deleuil et al. (2012) found that the uncertainty of the 
predicted Umb darkening coefficients was ±0.0168 when they 
propagated the stellar parameter uncertainties to the Umb dark- 
ening prediction. In most of the cases the stellar paramater un- 
certainties cause less than +0.03 uncertainty in the prediction of 
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Fig. 1. A comparison of 5 (3) different limb darkening models. 
We plotted only the coefficients of the quadratic limb darkening- 
law. From top to bottom one can see the models for different 
metallicities ([M/H] - -0.5, 0.0, 1.0). For the sake of cor- 
rect comparison, we used the same logg and turbulent velocity 
V, values for all the curves. Solid lines represent the ui coeffi- 
cients, dotted lines the U2 coefficients. The model of Sing (2010) 
is denoted by cyan lines. Red lines show the models of Claret 
& Bloeman (2011), which are based on the ATLAS synthetic 
spectra and their 'flux conservation method' (FCM). Green lines 
were obtained by them using the same spectra, but applying their 
'least square method' (LSM) for calculating the limb darkening 
coefficients. For solar metallicity ([M/H]=0.0), they also calcu- 
lated the limb darkening coefficients with these two methods, but 
from the PHOENIX spectra, too (violet and magenta lines). The 
figure only hows the CoRoT white light passband coefficients; 
the discrepancies are similar for Kepler's response function and 
for other photometric systems. 



the Umb darkening coefficients, (iii) Although the spectroscopic 
effective temperature determination is not significantly affected 
by stellar spots in the case of a solar-like activity, the stellar sur- 
face effective temperature determination can be systematically 
affected when the spot coverage is about 10-20% (A. Hatzes, 
priv. comm., Ribas et al. 2008), and this can lead to systematic 
errors in the limb darkening predictions and to underestimation 
of the uncertainties of such a prediction. Thus, the uncertainties 
mentioned in point (ii) might be underestimated. It is quite diffi- 
cult to characterize these kinds of error sources, but fortunately 
they have an impact only for very active stars, which are usually 
avoided by planet hunters. (A remarkable exception is CoRoT- 
2.) (iv) The inhomogeneties of the stellar surface temperature 
distribution (e.g. spots, faculae, gravity darkening-effects, etc.) 
are not included in the prediction of limb darkening coefficients. 

Now we consider the uncertainty in the planet-to-stellar ra- 
dius ratio caused by the uncertainty of the limb darkening coeffi- 
cients. After straightforward calculation we have (cf. Eq. (Bl)): 



AA: 
T 



1 / Auid + Au26^ 2Ami + Am2 

2 \ 1 - ui6 — U26^ 6 - 2u\ — U2 



(1) 



where we abbreviated 6-1 -fi, and Ax means the uncertainty of 
quantity x. The variation in the relative error in the radius ratio 
of the star and the planet from Eq. (1) is shown in Figure 2. 

First we consider the temperature region 5000K < T < 
7500 K. Here the different limb darkening theories agree more 
closely. Substituting the uncertainties of the theory discussed in 
the previous section, we can easily conclude that the uncertain- 
ties in the relative radius ratio of the planet and the star varies be- 
tween 1% and 10%, but it can be around 8% at the high temper- 
ature. Below 5000 K, where the curves of different theories are 
very divergent, we have uncertainties of up to 20% in the radius 
ratio. These uncertainties only come from the inconsistencies of 
the theoretical limb darkening calculations. The aforementioned 
two additional factors increase these uncertainties further, but 
one needs a case-to-case study to characterize them. 

These uncertainties are unacceptable if radii and therefore 
internal structure of the planets are being studied, because to 
distinguish between different planet models and to study them 
in detail we require +1% precision in planetary radii below 
five Earth-masses. The consistency of the limb darkening-tables 
should be improved in all temperature regions and should be 
checked via observations, to obtain more concordant limb dark- 
ening tables in the near future. 

Since we do not have reliable limb darkening theories it 
seems questionable to fix the limb darkening coefficients for 
transit light curve analysis. For instance, CoRoT and Kepler are 
able to determine the transit depth with a precision of 10 in 
general or sometimes even better, which can be translated to 

0. 5% relative error for the radius ratio ^ - in the absence of limb 
darkening. Since the error due to the poor knowledge of the limb 
darkening may be as big as 20% at low host-star surface temper- 
atures (r < 5000K), this can cause ~ 40 - 50 times larger error 
in the radius than the error stemming from the quality of pho- 
tometry alone! This error source is also two to ten times bigger 
at the higher host-star's surface temperatures {T > 5000/^) than 
the error stemming from the photometry. This is the case when 
one fixes the limb darkening coefficients during the transit light 
curve fit procedure. 

We add that these numbers come from the optimistic case; 

1. e., uncertainty in the stellar radius can be neglected. But this 
is not the case. The aforementioned errors are the errors stem- 
ming from the planet-to-stellar radius ratio determination alone. 



4 



Csizmadia et al.: Transit parameters, stellar spots and Umb darkening 





3000 4000 5000 6000 7000 8000 3000 4000 5000 6000 7000 8000 



Fig. 2. a): The difference of the theoretically predicted limb darkening coefficient ui as a function of the stellar eflFective surface 
temperature between Sing (2010) (S) and different models of Claret & Bloemen (201 1). F means the flux conservation method, L 
means the least square method of Claret & Bloemen (201 1). A and P denote the ATLAS and the PHOENIX models that they used 
for creating of their limb darkening tables, b): the same for the U2 quadratic coefficients, c): the absolute value of the maximum 
differences at a given temperature between Sing (2010) and some of the Claret & Bloemen (201 1) models. The curves in this panel 
are upper limits, because we calculated the differences between the different limb darkening tables in a conservative way; i.e., we 
searched for the biggest difference between Sing's table and any of the Claret & Bloemen's table at the given temperature, d): 
visualization of the meaning of Eq. (1): relative error of the radius ratio of the star and the planet k = RpUmetlRstar vs effective 
surface temperature of the stellar surface when the limb darkening coefficients are fixed. The Ami, Am2 uncertainties are taken from 
panel c. The two curves are valid for a rather central (b - 0.2) and for a rather grazing (b - 0.8) transit. For these plots we used 
the CoRoT passband limb darkening coefficients of Sing (2010) and Claret & Bloemen (201 1), and we assumed logg - 4.50, solar 
metallicity and a turbulent velocity of 2km/s. However, the main result does not change significantly for other input parameters or 
passbands. The uncertainty of the parameters were estimated as the difference between the two Umb darkening tables. In certain 
temperature regions we are far from the required precision in the radius ratio determination. 



When this radius ratio is transformed to absolute dimensions of 
the planet via RpUmet = kRstar, then the error bars should be in- 
creased more accordingly, by the uncertainties in the stellar ra- 
dius. The stellar radius is usually obtained by comparing stellar 
evolutionary models to the spectroscopically measured quanti- 
ties, such as Teff and logg of the star. It seems that the most 
careful spectroscopic studies of the host stars are able to yield the 
stellar radius with a precision not better than ~ 5% (e.g. Bruntt et 
al. 2010; see also the discussion in Torres et al. 2012). However, 
the final error budget also depends on the uncertainties and sys- 
tematics of the stellar evolutionary models used. For instance, 
the radii of stars below 1 Mq, predicted by current evolutionary 
models, are not supported by observations (see e.g. Clausen et 
al. 2009 and references therein). It is beyond the scope of this 
paper to characterize these kinds of error sources, we hmit our 
study only to the effect of Umb darkening on the planet- to- stellar 
radius ratio, and therefore we refer to Torres et al. (2010), where 
the uncertainties of the used stellar models and the derived stellar 



parameters are discussed. We also mention that future projects, 
like PLATO, are expected to provide the stellar masses and radii 
with less than couple of percent relative error via asteroseismo- 
logical studies. Details of this latter method and its uncertainties 
are discussed e.g. in Catala (2009). 

4. Adjusting the limb darlcening coefficients 

Here we study how precisely our fitting procedure is able to de- 
termine the different transit parameters. For this purpose we cre- 
ated 2000 synthetic light curves with the subroutines of Mandel 
& Agol (2002). We assumed a circular orbit for all of these 
curves. The number of the simulated photometric data points in- 
side the transit were selected randomly between 500 and 5000. 
Then we added Gaussian-type random noise to these Ught curves 
with zero mean and Icr scatter, and we chose randomly cr such 
that our final light curves have S/N between 1 and 1000. The 
S/N is defined as follows: S /N = k^jcr. We did not study the ef- 
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feet of red noise, systematic errors, etc., since our aim is to find 
what minimum accuracy is needed for precise parameter deter- 
mination. The input parameters of these synthetic light curves 
were randomly chosen between the following limits: b - 0...1, 
k = 0...0.3, a/R, = 5...50, m+ = 0...1 and u- = -0.1... 1.0, re- 
spectively. Here, a is the semi-major axis, Rs the stellar radius - 
a/Rs are often called scaled semi-major axes -, b is the impact 
parameter, and m+ = mi + U2 and m_ = mi - U2. The use of these 
combinations of the limb darkening coefficients was suggested 
by Brown et al. (2001) and they aim to avoid degeneracies be- 
tween Ml and U2- We selected u+ and u- randomly between the 
aforementioned limits. This means that they are not related in 
this test to any stellar properties or theoretical considerations, 
and thus there is no link between them; i.e., they are indepen- 
dent of each other in this test. 

In the next step we modelled these artificial light curves with 
our fitting procedure, which is based on a genetic algorithm 
(Geem et al. 2011). The free parameters were the impact pa- 
rameter b, k, a/Rs ratio and m+. We fixed the value of m_ at 13 
difi'erent values at -0.3, -0.2, -0.1, 0.9, 1.0, and then we se- 
lected the best solution (according to the ^-values). Then, this 
best solution was refined by allowing m_ to vary in a narrower 
range (±0.1). After that we plotted the standard devations of the 
difl^erences between the value obtained from the light curve so- 
lution and the original input value, normalized to the input value 
vs the S/N ratio as it is plotted in Figure 3: 

/ Imodelled value - input valuel \ 

s = STD DEV : . (2) 

\ mput value / 

This s-value has been plotted in the different panels of Figure 3 
against the S/N of the simulated light curves. We draw conclu- 
sions from this figure as follows. We need at least S/N=25 and 
6 to determine the hmb darkening coefficients and the scaled 
semi-major axis with ±5% uncertainty. To determine k with 1% 
accuracy, we need at least S/N=~50. To determine b with 10% 
relative accuracy, we need S/N=~ 100 at least. We also found 
that when b > 0.85 the solution becomes unstable and far from 
the input values, so we removed them from this analysis. 

5. Effect of stellar spots on the limb darkening 

The result of a fit with adjusted limb darkening coefficients can 
lead to unreasonable limb darkening values, e.g. Southworth 
(201 1) accepted a worse solution for CoRoT-13b in the sense of 
the quality of the fit (measured by the;^^-values), because when 
he adjusted both hmb darkening coefficients, he obtained 'un- 
physical' values, i.e. far from the theoretically predicted range, 
sometimes causing 'limb brightening' instead of limb darken- 
ing or negative fluxes at the edge. What can we say about such 
cases? 

We would like to point out that we employed an indirect as- 
sumption in the previous sections: the effective temperature of 
the stellar surface is assumed to be the same at each point of 
the surface and therefore the stellar surface brightness distribu- 
tion is only modified by limb darkening. The same assumption 
is used when the theoretical limb darkening coefficients are cal- 
culated, and almost every transit Ught curve investigator uses the 
same assumption. However, there are several factors that modify 
the local effective temperature on the stellar surface. Such ef- 
fects are e.g. stellar rotation (i.e. gravity darkening, see e.g. von 
Zeipel 1924; Twigg & Raffert 1980; Bames 2009, Espinoza & 
Rieutord 2011; Claret 2012), reflection effect (although we can 
neglect this for the star in a star-planet system), and stellar ac- 
tivity phenomena: spots, plage-area, faculae and flares. Some of 



these are transient phenomena, e.g. the flares, but others have 
similar or much longer timescales than the length of the transits 
(e.g. the spots and faculae). They can be regarded as an addi- 
tional, but constant surface brightness distortion during a transit. 
Constant here means that we assume that the number, tempera- 
ture, size, astrographic longitude, and latitude of the distortion(s) 
do not change around and during a transit. The constancy of the 
aforementioned quantities are reasonable. 

These distortions mean that the local effective surface tem- 
perature where they occur is lower for stellar spots and higher 
for faculae than that of the undistorted stellar surface. As is well 
known, the limb darkening coefficient is different for different 
surface temperatures. This means that the apparent stellar disc 
cannot be characterized by only one effective surface tempera- 
ture value, and that is why the surface brightness distribution 
cannot be characterized by only using the limb darkening coef- 
ficients due to only one Teg. 

Djurasevic (1992) already took into account that the undis- 
torted and the distorted areas have different limb darkening co- 
efficients, when he developed a code to analyse eclipsing bi- 
nary star light curves. The corresponding limb darkening coef- 
ficients were taken from theoretical approximations available at 
that time. We would need a similar approach for the transit mod- 
elUng. However, we suggest avoiding exactly the same approach 
for the following reasons. 

a) Different tables give different Umb darkening coefficients, 
which is why the results can be different, and maybe the 
derived spot parameters will depend on the assumed limb 
darkening coefficients; 

b) The spot distribution is not weU known, and the spot solu- 
tion is often highly degenerated so the reliability of the spot 
solution can be criticized if only one-colour photometry is 
available. Consequently it is hard to include the spots in a 
modelUng tool. 

5.1. Effect of a single active area 

To investigate the effect of stellar spots and faculae (hereafter we 
call both spots, including spots that can have higher temperatures 
than the surface, i.e. "bright spots", describing the characteristics 
of faculae that they have a higher temperature than the normal 
photospheric area), we classify the spots as follows. Spots that 
are not eclipsed by the planet we define as Type I spots and spots 
that are eclipsed by the planet as Type 11 spots. 

The effect of Type 11 spots can be taken into account because 
the size and location of the spots can be determined with high 
accuracy, and that is why removing of their contribution from 
the light curve causes no problem, at least in principle (Silva- 
Valio & Lanza 201 1; Sanchis-Ojeda & Winn 201 1; Desert et al. 
2011). 

Type 1 spots are more complicated. Those that cause light 
curve modulations can be removed safely with baseline fits and 
baseline corrections using the out-of-transit points in the small 
vicinity of transits. But there are also spots that cause no light 
curve modulations, e.g. a polar spot. In addition, Jackson & 
Jeffries (2012) propose, based on spectral evidence, that a con- 
siderable amount of stars exhibit numerous (up to ~ 5000), 
small (not bigger than 2-3 astrographic degrees) dark spots that 
are axysynmietrically distributed on the steUar surface, reaching 
~ 50% spot-coverage, which led to no observable light curve 
modulation over the current detection limit in spite of their high 
activity level. If this is true then we can expect that even pho- 
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tometrically quiet stars can produce strange limb darkening co- 
effients as we will show. 

Figure 4 illustrates the situation. When the planet moves 
in front of an unmaculated star, the observable maximum flux 
change will be 



AF _ j^ l- Mlo(l - IJp) - M20(l - flpf 

— -r 1 



£10 _ jf20 

3 6 



(3) 



where AF is the flux decrease due to the planet transit, Fq is the 
total emitted flux of the star in our direction, uio and M20 are the 
limb darkening coefficients of the unmaculated steUar surface 
and index p stands for the planet. 

If we have one circular Type I spot on the stellar surface, then 
this can be modelled by replacing the stellar flux with the spot's 
flux at the place of the spot. Then, we use the new Fstar value 
instead of Fq in Eq. (3). In this case, the total flux of a spherical 
star with only one Type I spot can be written as (if the spot is 
small enough) 



nR^f^^kil - Mio/3 - M20/6) 



(4) 



Index s stands for the spotted area. The first term on the right 
hand side describes the total flux stemming from the stellar disc, 
while the second term describes the effect of the spot: we remove 
a certain small part of the stellar flux and replace it by the spot's 
flux. The argument of the Lo function is which we can calcu- 
late easily from the spot's position. This equation also shows that 
in the case of a dark spot, Fs,ar is smaller and the transit depth 
becomes larger. For faculae, which are brighter than the surface, 
Fstar become larger and hence the transit depth will be smaller. 

When we model the transit, we replace this real star with a 
hypothetical star of an unmaculated surface. Although the tran- 
sit depth will change due to the presence of stellar spots, it is 
easy to show that the radius ratio k remains unchanged and pre- 
cisely determinable if one chooses the effective Umb darkening 
coefficients appropriately. 

We use the following abbreviation (cf. Eq. (4)): 

N 2 

^jll^RsjUoLoiUm, "20. A'sj) - hjLD{U\sj, U2sj,IJ.sj)] 

^ = — n 

where j stands for the jth spot. Then we have a step-by-step 
calculation of the observable light loss during the transit. 
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+ a, 



a = 



3(1^-6(1^) 



(7) 



This a is constant during the transit, because it does not depend 
on the planet-star distance and we assumed above that the spot is 
fixed in space, i.e. the star rotational period is much longer than 
the transit length. 



Finally the basic equation which relates the observed light 
loss during a transit, the planet-to-stellar radius ratio, stellar spot 
effects and limb darkening to each other is 



AF- 



1 



"le/f 
3 



6 



(8) 



Comparing the equations of this section to each other we can see 
that the newly defined effective limb darkening coefficients are 
related to the values vaUd on the unspotted stellar surface such 
that: 



Uuff - 



1 -A' 



U2eff - 



U20 
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(9) 



From this mathematical analysis of transit light curves it is obvi- 
ous now that spots acts Uke a contamination source in the aper- 
ture and this is valid during the entire transit. If an unresolved 
star or other light source (also called "contamination source") 
contributes to the observed flux, then that decreases the observed 
transit depth. This extra light source is routinely removed from 
the light curve in general. Certain spots acts as a negative a. In 
total, this kind of contamination-like source can be either posi- 
tive or negative, but in the transit fits we did not include negative 
contamination. 

It has been already recognized, that Type I spots can act as 
contamination sources, but our result implies also that limb dark- 
ening coefficients can also change which is a new result. In ad- 
dition to the arbitrary sign of a, these latter two equations jus- 
tify that contamination should be an additional free parameter in 
the transit fit which would enable to remove all the spot-effects 
which may remain after the baseline-correction since they might 
not cause light curve modulation in and out of transit. It also 
shows that no detailed spot modelUng is needed to fully under- 
stand the transit light curves if no spot-crossing occurs, even if 
polar or any other spots are present. The problem of this case can 
be solved in the following simple and robust way: limb darken- 
ing coefficients and contamination factor should be free param- 
eters. Last but not least it shows that the observed limb dark- 
ening coefficients can be highly different from the theoretically 
predicted values in the presence of spots, even if they do not pro- 
duce light curve modulations. The uniqueness of the transit light 
curve solutions are provided by the fact that the length of the 
ingress and egress phases, as well as the full transit length, are 
also related to the size ratio of the star and the planet (see e.g. 
Winn 2010). 

This analysis is only valid for small spots. The analysis can, 
however, be generalized easily for numerous and larger spots, as 
we did with the summation above. It is also easy to extend the 
analysis to bright spots (i.e. faculae). If one takes more terms in 
the expression of the limb darkening laws into account, then the 
effective Umb darkening coefficients will be transformed in the 
same way as the first two terms, and it is then clear how one has 
to extend the a-function with additional limb darkening terms. 

We also note that for this spotted case the same applies as for 
the previous unspotted case in Section 4: if the limb darkening 
coefficients are fitted, then we can determine their exact values if 
they fulfil the S/N requirements. We have to keep in mind, how- 
ever, that we measure the effective limb darkening coefficients 
in this case and not the ones related to an unmaculated stellar 
surface. 

Recently, combined optical and near infrared observations of 
transits were proposed by Ballerini et al. (2012) to minimize the 
effect of spots on the transit parameters. It was proven therein 
that it is possible to derive precise planetary radii from this kind 
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of multicolour observation, if the temperature of spots and the 
theoretical Umb darkening coefficients are known. Here we have 
showed that the same result can be reached without spot parame- 
ters and without the knowledge of the theoretical limb darkening 
coefficients (which are practically unkown today). Our method 
is not only simpler, but it does not require multicolour photom- 
etry. Both methods have the advantage of being able to provide 
the exact and precise planetary parameters in highly active cases 
when the spots are not crossed by the planet, and therefore we 
have limited information about the stellar surface. 

5.2. Effect of more than one active area 

The direct observations of the Sun's surface, the Doppler- 
imaging of other active stars and their photometric spot studies 
show concordantly that, in general, more than one active area is 
present on the stellar surface at the same time. The results of the 
previous section can be generalized easily by simple summing 
of the effects of several individual active areas (cf. Eq. (5)). 

Studies of the Sun show that the total area of the faculae 
depends on the area of the dark spots. According to Chapman et 
al. (1997), the following expression is a good approximation for 
the Sun 

'^faculae ~ ^^^spot, (10) 

but the conversion factor of the total area ratio shows small sec- 
ular variations during one solar cyle (Chapman et al. 1997). It is 
reasonable to assume that other stars have similar relationships 
between their total area ratios, although the conversion factor in 
Eq. (10) can vary. For instance. Chapman et al. (201 1) find that 
the conversion factor varies from cycle to cycle, and for solar 
Cycle 23 it was 42 instead of the 17 that is vaUd for Cycle 22. 
More studies are required to determine this factor and its depen- 
dence on the stellar properties. Foukal (1998) shows that there 
is an anticorrelation between Agpo, and A faculae in the case of the 
Sun and he argues that this can be extrapolated to other stars. 

The total spot coverage varies over a wide range. The Sun 
has a spot coverage of 0.1% only at maximum. Several observa- 
tions suggest that many of the normal stars can have 1-2% spot 
coverage, but there are quite a few stars showing 10-60% spot 
coverage on its surface (Schrijver & Zwaan 2004 and references 
therein; this coverage refers to the total area of the spots rela- 
tive to the whole stellar surface). The most active stars seem to 
be located in close binary stars and giants, and they are not so 
important from the point of view of the presently known transit- 
ing exoplanets. But some exoplanet hosting stars (like CoRoT-2, 
CoRoT-6, CoRoT-7 etc.) seem to be very active stars. E.g. Lanza 
et al. (2009) established ~ 7% spot coverage for CoRoT-2. 

The changes in the limb darkening coefficients have dramatic 
conclusions for the errors in planet-to-stellar radius ratio if one 
fixes the limb darkening coefficients according to tabulated val- 
ues and neglects the effect of spots on the limb darkening. If one 
accepts the conclusions of Jeffries & Jackson (2012), namely 
that ~ 50% of the surface can be covered by small spots in 
an active star without observable photometric modulation (but 
spectroscopic activity indicators still imply an active nature of 
this hypothetical star), then the A-function defined by Eq. (5) 
can have a value of ~ 0.36. This causes via Eq. (9) the observ- 
able limb darkening coefficients to increase by a factor of 1.56 
- if we do not take into account the role of faculae, which may 
be negligible according to the extrapolation suggested by Foukal 
(1998). If the host star is less active, we take 5% spot coverage 
as an example, the A-function will cause a change of ~ 0.037 in 



the limb darkening coefficients, which will lead to 9% error in 
the radius ratio (cf. Eq. 1)). If the star has only 0.5% spot cov- 
erage, then the A-function will cause a change of only 0.004 in 
the limb darkening coefficients, but this still causes 1 % error in 
the radius ratio. Thus, 0.5% spot coverage can be considered as 
a limit: below it one can neglect the change in the limb dark- 
ening coefficients owing to Type I stellar spots. These changes 
in the observed effective and the tabulated limb darkening coef- 
ficients are plotted in Figure 5. (For this estimation we used a 
rather grazing eclipse. For central eclipses the effect is less but 
still considerable.) 

We note that when a host star of a transiting planet is highly 
active, like CoRoT-2, then it is not enough to take an appropri- 
ate contamination factor into account to change the transit depth 
or to carefully model the crossing of Type II spots, including a 
baseline correction to remove all the contributions of spots to the 
photometry. It is also necessary to keep the limb darkening co- 
efficients as free parameters or to modify them according to Eq. 
(9), because the modification of the limb darkening coefficients 
will also cause non-negUgible effect on the planetary radius. 

Thus we conclude that stellar activity modifies both the ob- 
servable effective limb darkening coefficients, as well as the ob- 
servable transit depth. Then, if spots and faculae are present, 
one cannot fix the limb darkening coefficients corresponding to 
the stellar effective surface temperature. It is preferable to leave 
them as adjustable parameters, because it is difficult to predict 
the effect of the spots on the limb darkening coefficients. We also 
established that the effective limb darkening coefficients (note 
that almost all transit Ught curve modellers use these) can have 
unphysical values because of spots and faculae. Of course, 'un- 
physical' is only the correct word if, and only if, we do not com- 
pare the spotted star to an unspotted model calculation. It is also 
clear that in such a case, we also can determine with good ac- 
curacy (< 1%) the radius ratio of the star and the planet, the 
impact parameter, and the aIRs ratio and the limb darkening co- 
efficients, but we just keep in mind, that in the case of several 
spots and faculae, the limb darkening coefficients yielded by the 
modelling do not correspond to the hmb darkening coefficients 
of an urmiaculated star with the same effective temperature, but 
instead correspond to the effective limb darkening coefficients. 

6. Determinability of individual limb darkening 
coefficients 

Here we study the determinability of the specific limb darkening 
coefficients. Since many combinations of Umb darkening coef- 
ficients are able to reproduce the same surface brightness distri- 
bution within the errors of the observations, it seems likely that 
the light curve fitting procedures reproduce the intensity distri- 
bution and not the true values of the coefficients. In other words, 
the inversions of the polynomial expressions of the limb darken- 
ing laws incorporate many degeneracies. That is why we are able 
to determine all the transit parameters with good accuracy and to 
reproduce the corresponding stellar surface brightness distribu- 
tion within the errors of the observations, but we are not able 
to determine the individual values of the Umb darkening coeffi- 
cients. From the point of view of exoplanet parameter determi- 
nation this does not present any difficulties. To illustrate this, we 
rewrite Eq. (1) for the following analysis as 



The current aim is to know the planetary radii with 5% precision 
(Valencia et al. 2007), although the future goal is 1% (Wagner et 
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al. 201 1). To determine k, say, with 5% precision, it is enough to 
know the limb darkening with a precision of 10%, which follows 
from the Unear error analysis 



^k _ 1 l^F ALd 



(12) 



because the photometric error is practically negligible for space- 
based photometry beside the error steming from limb darken- 
ing's errors (cf . Section 3). This means that the transit light curve 
modelling codes can produce quite different limb darkening co- 
efficients than what we expect. Now we aim to determine these 
differences. 

We denote the expected limb darkening coefficients by m,, 
and the difference between the measured and the expected val- 
ues by ySj. The relative deviation between measured and expected 
stellar intensity profiles is denoted by and the maximum al- 
lowed relative deviation - occurs at the still unknown 6 - is de- 
noted by Cmax, so we can write 



Ld{U\,U2, ...,S) - Ld{Ui +p\,U2 +P2,-,S) = emax- 



(13) 



Although the analysis can easily be extended to higher orders, 
we cut it at second order because this is generally used for tran- 
sit analysis (but as can be seen, the conclusion does not change 
significantly if one extends ffie calculation to higher orders). 

We chose a solar-like star to illustrate ffie operation of the 
j6,-numbers. We plotted the deviation |e| for different cases in 
Figure 6. In ffiat figure, ffiere is a rather wide range for yS, pa- 
rameters ffiat produce the same hmb darkening function within 
+e tolerance range. It is easy to construct an Lo function, for 
example, of limb darkening coefficients of u\ + Pi - 0.8 and 
U2+ Pi = -0.1 and this produces the same results wiffiin the 
tolerance range as ffie correct values ui = 0.41 and U2 = 0.26. 
This example is presented in Figure 6. Along the sides of the 
paralellograms one can easily construct the same quality of ffie 
light curve fit (cf. Eq. (13). In Figure 7 we plotted anoffier ex- 
ample: ffie stellar surface brightaess distribution for different j8, 
numbers. As one can see ffiere is a high tolerance range, so one 
can choose very different limb darkening combinations from ffie 
theoretical predictions that produce a similar surface brightness 
distribution wiffiin ffie tolerance range. This example shows that 
we are able to produce acceptable D(r/R - /u) functions within 
the tolerance ranges wiffi quite different limb darkening coeffi- 
cients. If the tolerance range changes due to diflFerent precision 
requirements, the acceptable range of ffie limb darkening coeffi- 
cients will also change. 

One can envisage ffiat a transit Ught curve modelling code 
may take different numbers as input for numerical reasons or 
owing to data points/error distributions than the theoretically ex- 
pected values. This can be also the case if stellar spots modify 
ffie hmb darkening coefficients. That is why ffie absolute values 
of ffie limb darkening coefficients are not important but the shape 
of the radial intensity distribution that they define is (see Figure 
7 for such a comparison). 

This analysis gives hints and arguments that what we de- 
rive from ffie transit hght curve are not the true individual values 
of ffie hmb darkening coefficients but some raffier complicated 
combination of them. The limits of the potentially wide ranges 
depend on ffie actual impact parameter and photometric preci- 
sion (i.e. ffie range of the tolerance factor e), so that to give the 
hmits requires case studies. 

From Figure 6 it is also clear ffiat the fitted limb darken- 
ing coefficients are also correlated with each other. This be- 
haviour has already been observed by numerical experiences 



in Brown et al. (2001) and in e.g. Deleuil et al. (2008) and 
Csizmadia et al. (201 1). Brown et al. (2001) and Pal (2008) sug- 
gest fitting the combination u+ = u\ + U2, u- = u\ - U2 and 
M+ = Micos40° + M2sin40°, m_ - micos40° - M2sin40° in- 
stead of fit Ml and M2 themselves because ffiey are less correlated 
to each other as previously assumed. As is clear from Figure 6 
this rotational-hke transformation can help better constrain one 
of ffiese combinations, but the offier will still feature a value in a 
wide range and the problem cannot be surpassed in this way. The 
values that look 'unphysicaF at first sight may describe the radial 
intensity distribution as well as do other values. We emphasize 
again ffiat the absolute values of the fitted limb darkening co- 
efficients are themselves not important but ffie radial intensity 
profiles are. These profiles, taking the tolerable deviations into 
account, can be produced with a very wide ranges of limb dark- 
ening coefficients. That is why the basis of the analysis of the 
results are not comparison of the resulting limb darkening co- 
efficient values themselves to some theoretical table values, but 
the comparison of the resulting radial intensity profiles of the 
predicted wiffi the observed ones, a diagram similar to Figure 7. 

7. Conclusions 

In this paper we have attempted to understand ffie impact of limb 

darkening on ffie planetary radius, to present ideas as to which 
effects can modify the theoretically predicted limb darkening co- 
efficients, and to asnwer why some authors have highly different 
observed hmb darkening coefficients from the theoretical pre- 
dictions. 

We concluded that concurrent limb darkening tables do not 
yield consistent results in several temperature- and metallicity 
regions, and we showed that the inconsistencies will lead to se- 
rious defects in ffie determination of ffie planet-to-stellar radius 
ratios if one fixes ffie limb darkening coefficients according to 
ffiese tables. We also proved with numerical experiments ffiat 
one can safely fit ffie limb darkening coefficients. 

We investigated how the presence of stellar spots - both dark 
spots and bright faculae - modify the observable transit depffi 
and ffie limb darkening coefficients. We found ffiat below 0.5% 
spot coverage, the change in the effective limb darkening co- 
efficients causes no significant error in the measured planet-to- 
stellar radius ratio, but over ffiat coverage ffiis cannot be ne- 
glected, and neglecting ffie change in ffie hmb darkening coeffi- 
cients ffiemselves has a significant change in ffie radius ratio. 

We pointed out that the radial intensity profiles of the stel- 
lar discs are not very sensitive to the specific values of ffie limb 
darkening coefficients, and ffiey are anticorrelated to each offier 
in the case of the quadratic limb darkening law (cf. Figure 6), 
and a similar statement can be made about higher order limb 
darkening laws. Consequently, the observed values of the limb 
darkening coefficients can be very different from the ffieoreti- 
cally tabulated ones even if ffiey produce ffie same radial inten- 
sity profiles within the tolerance. 

Limb darkening is still a problem for transit modelling. On 
ffie one hand, limb darkening is a challange for stellar astro- 
physics. On the offier hand, we can replace the hmb darken- 
ing coefficients of an unmaculated star wiffi ffie appropriately 
weighted effective limb darkening coefficients for the modelling 
Our final conclusion is that the transit parameters can be deter- 
mined even if our knowledge of limb darkening and stellar spots 
is not satisfactory. Our work will make it clear that one has to 
fit ffie limb darkening coefficients and it is not self-evident that 
the resulting limb darkening coefficients should agree with the 
ffieoretical predictions. Indeed, ffiese can be very different. Our 
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results suggest not only that fitting the limb darkening coeffi- 
cients will produce more reliable solutions but also that their ad- 
justment will not affect the accuracy of determining other tran- 
sit parameters. However, predicting the effective limb darkening 
coefficients is difficult, because it requires a reliable and detailed 
spot solution. Moreover, these can vary from transit to transit as 
the spots evolve. If we fit them, then we can determine their val- 
ues, too, but then such effective limb darkening coefficients have 
no direct link to the theoretical limb darkening tables. When a 
statistically significant sample of high-quality, multicolour pho- 
tometric analyses of transits become available, then the compar- 
ison of the resulting limb darkening coefficient to the available 
tables will become useful to characterize the discrepancies be- 
tween theory and observations better This can contribute signif- 
icantly to improving our knowledge of stellar atmospheres. It is 
also important to simultaneously determine the activity level of 
the star and the limb darkening coefficients 

It is also clear that further detailed studies are required to 
completely understand the limb darkening phenomenon and its 
interaction with e.g. gravity darkening. Such studies are very 
useful for planning such projects as PLATO or for understanding 
the present-day ground- and space-based measurements more 
deeply. 
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Appendix A: The limb darlcening laws 

Limb darkening is a stellar atmosphere phenomenon that is re- 
viewed in detail in several textbooks (see e.g. Kallrath & Milone 
2009). The most general formula to describe the effect of limb 
darkening is 

I(y) = loLoiui, U2, Ui, ...,n) = (A.l) 
7o(l-S;^iM/l-cosr)^), 

where Uj are the so-called limb darkening coefficients, y the an- 
gle between the normal vector of the surface point and the di- 
rection to the observer, /q and I(y) are the surface brightnesses 
(or intensities) at the appearant centre of the stellar disc and in a 
certain point of the stellar surface characterized by the angle y, 
and Lo is the limb darkening function defined by Eq. (A.l). We 
also use the common notation jj. - cos y. 

Many studies of eclipsing binaries have indicated that at least 
two terms in Eq. (A.l) should be kept to fit the light curves sat- 
isfactorily (Twigg & Raffert 1980; Alencar & Vaz 1999; Albrow 
et al. 2001; Claret 2000, 2008, Southworth et al. 2007). That is 
why, in the exoplanetary research field it is widely accepted to 
use the truncated formula of Eq. (A. 1) 

LD = l-'^%iUjii-cosyy, (A.2) 
which can easily be rewritten as 

Ld = I - iui+ U2) + (ui + 2u2)n ~ U2H^. (A.3) 

There are different approaches describing the limb darkening ef- 
fects, including logarithmic or square root terms of /z, see e.g. 
Claret (2000) or van Hamme (1993) for an overview. 
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Appendix B: Transit depth and limb darlcening 

It is straightforward to show in the small-planet approximation 
that the transit depth can be written as 



AF ^ j^ l - ui(l - fi) - U2(l - nf 
F 1 _ !fl _ !!2. ' 

^ A <c 



where 



and 



„ d2 /i "1 "-A 



(B.l) 



(B.2) 



(B.3) 



also hold, and it is trivial to extend these expressions to higher 
orders of the limb darkening approximations (cf. Kallrath & 
Milone 2009 and Mandel & Agol 2002). In these expressions 
AF is the light lost due to the transit; F is the stellar flux out of 
the transit part of the hght curve; k — R planet star 

is the ratio of 

the planet radius RpUmet and the star's radius Rstar- In addition, in 
the moment of the maximum light lost, ji - Vl - where b is 
the impact parameter. If we are interested in the light loss at any 
arbitrary time inside the transit, then yu = Vl - (5^ where S is the 
sky-projected distance measured in stellar units. (For details on 
how to determine it for any arbitrary time see Kallrath & Milone 
2009.) 

Taking the typical values of the limb darkening coefficients 
into account (see Section 2 or the references mentioned in 
Section 1), one can easily see that limb darkening will modify 
the true depth by up to a factor of ~ 2 relative to the case where 
one takes the transit depth to be proportional only to the radius 
ratio - therefore the eflsct of Umb darkening is not negUgible. 
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Fig, 3. Results of the light curve modelling test. On the x-axis 
the S/N quantity is defined in the text, y-axes are the relative un- 
certainties in radius ratio k, u+ combination of the two quadratic 
limb darkening coefficients used, in the scaled semi-major axis 
(for brevity we used the notation of a only instead of a/Rs in the 
third panel from top) and in the impact parameter b, from top 
to bottom, expressed in percentages. The solid lines show the 
relative uncertainties defined as the standard deviations of the 
difl'erences of the modelled and the input values normalized to 
the input values at certain S/N- values. 
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Fig. 4. Illustration of the effect of Type I spots. Left: the planet 
crosses an unmaculated star that is characterized with some Umb 
darkening coefficient uq. Right: the planet crosses the appareant 
stellar disc of a spotted star, where the spots and the planet have 
different impact parameters, as well as the stellar photosphere 
and the spots have different limb darkening coefficients (mq, m,). 
Grey area is the planet, black ellipses represent the spots. 




Spot coverage (%) 

Fig. 5. The x-axis is the total spotted area in percentage of the 
whole stellar surface area. The y-axis is the effective - i.e ob- 
served - limb darkening coefficients relative to the table value at 
the given stellar efffective temperature. For this figure we used 
Tsiar - 5775K and Tspo,s = 3775K, and the positions of the spots 
were chosen randomly on the visible hemisphere. The size of the 
spots were always the same, so higher spot coverage corresponds 
to larger number of spots. The Umb darkening coefficients were 
taken from Claret & Bloemen (2011) for /?-band, log^ = 4.5 
and at solar metallicity at these temperatures. Effects of faculae 
were not considered in this figure. 
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Fig. 6. The x-axis is the modified Unear limb darkening coeffi- 
cient xi = Ml -H ySi, the y-axis is the modified quadratic limb 
darkening coefficients = mi +y62- "i, "2 are fixed. The point 
in the middle shows the u\ = 0.41, U2 - 0.25 point, i.e. the 
original value of the limb darkening coefficients of the example 
presented in Section 8. Contours represent the quantity emax, i-e. 
the absolute value of the maximum deviation of the original Umb 
darkening function with mi,2 coefficients from the Umb darken- 
ing function with xi_2 (S run between and 1). 
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Fig. 7. The radial intensity profile Dir/R) as a function of the 

sky-projected distance r from the stellar centre (R is the stellar 
radius). As in all calculation in this paper, /i = -v/l - (r/R)^. 
The soUd line shows the effect of limb darkening D{r) - 
Lo{u\,U2,n)l{l - Ml/3 - M2/6) that - multiplied by - is di- 
rectly proportional to the light loss during a transit. The dashed 
lines show the tolerable ranges: between these lines the radial 
intensity distribution profile will produce a radius ratio k that is 
in the tolerance range of ±5%. The dotted line is an example 
of an acceptable radial intensity distribution profile with ui = 
0.82, M2 = -0.16, while the dot-dashed line is with u\ - 0.02, 
U2 = 0.6. 



